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ABSTRACT 


This thesis describes briefly the scope and nature 
of multi-echelon inventory systems where item repair is 
possible at the various levels. A model is created which 
describes many commonly encountered real world systems. 
feTenot NON -comparable and comparable costs are specified 
and a cost-effectiveness approach to the solution of the 
model is outlined. A nonlinear program is developed where 
the expected sum of comparable costs is minimized subject 
to a given level of fill rate (effectiveness) provided by 
the lower level stocking points to their customers. 

Beer eaters 26re tnen made to indicate possible soiu- 
Eton proceaures tor the program and characteristics and 


Mees) tOr Ghe COst—-errectiveness alternatives developed. 
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I. INTRODUCTION 


This thesis deals with a complex model of multi-~-echelon 
Maiventory Control in which the items are considered to be 
repairable. Whenever an item fails, the user demands a 
replacement and turns in a nonfunctioning unit called a 
carcass. An attempt is made to repair the carcass and 
[meurn tt to stock for future issue. Characteristics of 
items designated as repairable are: 1) such items are phy- 
Sically capable of repair, 2) the cost cf repair is less 
than procurement cost, and 3) initial procurement costs 
Pmemtrorge. Complicating this scheme is the fact that nearly 
el Ssuch models are of the multi-echelon type. That is, 
Several levels Of Gepair and stockage interact in compli- 
Cated ways which greatly increase the difficulty cof analy- 
Sis. Multi-echelon systems are characterized by a group of 
Stocking activities which receive customer demands and in 
turn demand resupply from a higher level. There may be any 
number of levels or echelons; however, in practice only two 
or three levels are usually encountered. As an example of 
a multi-echelon inventory system consider a large retail 
Store chain which procures and stores certain types of goods 
nationally and, when required, supplies regional warehouses 
Vile Weebimeicn SUSDlyY Jocal retail outlets. The objectives 
of the study of multi-echelon inventory control deal to a 


large extent with the development of those policies of 





erdering at all levels which produce some specified system 
performance, such as minimum expected system operating 
Costs. 

In addition to being complex multi-echelon repairable 
item inventory systems also often entail huge expenditures 
M=mmi@iey and Owaer resources. One of the largest multi- 
echelon inventory systems dealing with repairable items is 
managed by the U. S. Navy's Aviation Supply Office (ASO). 
Repairables represent 18.7% of the items controlled by ASO 
py type, £or a total of 66,000 separate items; but, almost 
60% by dollar value for an investment of neariy 3.1 billion 
Semlacs. The Aviation Supply Office spends 16 million 
dollars per month for repairs at the lower levels of repair 
alone, that is, simple or relatively inexpensive repairs, 
and 28 million dollars monthly for procurement of such 
items. Among ASO managed repairable items, 365,000 units 
Pee reworked in fiscal year 1971. Im view of these large 
investments in multi-echelon repairable item systems it is 
Met GQitticult to understand how a small improvement in the 
way in which these systems are managed can result in sub- 
Stantial savings. 

This thesis describes a model which approximates the 
Mee aeworking Structures of many repairable item, multi- 
echelon inventory systems. A set of costs and measures of 
erfectiveness is defined and its application to the model is 
discussed. The model is solved for minimum expected total 


System cost ordering policy subject to a speci ied level of 





effectiveness or performance. The solution is presented in 
the form of a nonlinear programming problem. Several 
remarks are then made in an attempt to outline possible 
solution procedures for the program and applications of the 


model. 





El. @2eEeMODEL 


The lower level of the multi-echelon repairable item 
inventory system considered in this thesis will consist of 
two bases which experience demands resulting from failures 
of items in use. When a demand is made on base i, a 
Carcass or not-ready-for-issue item is turned in at that 
base's repair facility. In practice it is the case that 
some of the carcasses which are turned in at each base 
Cannot be repaired at the base. It will be agsumed in this 
thesis that an item will be repaired at base i with proba- 
baa ty Pj- If repair at base i is possible the repaired 
item is returned to stock at base i after some time, rt... 
later. It is assumed that rt. is a random variable with 


peme general Gistribution and has mean rt Mies ee ais 


i - 
policy at base i may then best be described as an infinite 
server queue. If an item is deemed to be unrepairable at a 
base it is forwarded to a higher level repair facility, 
depot repair. Define the quantity of items on hand in 
Stock and in repair plus items on order less items ordered 
by customers but not supplied at base i to be the base i 
inventory position. When the base i inventory position 
fegemeoead level rf, a request is made to the depot for an 
Mi@mememel Unit which arrives a time t; later where t, 


Lona vaio Oof a cram@om variable. Since the base i inventory 


position only decreases whenever a demand is received and 





the accompanying carcass 2s not repatrable, an order is 
placed only when such a demand occurs. 

A single depot receives demands and carcasses from the 
bases when repair is not possible at base level. It is 
clear that some items will eventually be lost since their 
Geme7tion will have deteriorated beyond repair. For this 
reason it is assumed that carcasses received from base i 
may be repaired with probability dy ane If repair is met 
possible they are discarded. Items accepted for repair are 
assumed to be returned to depot stock after a constant 
repair time, rt. Define the guantity at depot both on hand 
and in repair to be the depot net inventory. Due to the 
attrition caused by items being discarded, the depot net 
Inventory position will fluctuate. Define the quantity on 
hand in stock and in repair plus items placed on order by 
ie Oecpec 1eeegenmcers placed by the bases to the depot but 
Pemailning UnLwEtea to be the depot inventory position. 

[een Che depot inventory position reaches a level fr, an 
Sieger iS “initiated by the depot for a quantity Q froma 
Pounce Of manufacture or higher level of supply which 
arrives after a constant lead time t. This policy is simi- 
lar to the well known (r,Q) model of single stock points 
which has many practical advantages as well as intuitively 
appealing logic and simplicity. Among its advantages is 

the ease with which it may be understood and applied to 
nearly any inventory requirement. It is necessary to con- 


Sider the items in depot and base repairs wher developing 





an inventory ordering policy. In this thesis this is 
accomplished by imcluding the quantity in repair in the 
inventory position. In this manner it is possible to 
insure that the quantity of items in the entire system 
remains bounded. The movement of items within the model is 
PeeoeGerea in Pig. I. A listing of parameters and vari- 
ables used in the description and solution of the model is 


included in Table I. 


A. MAJOR ASSUMPTIONS 
Major assumptions about the model required in its 
solution are: 


base 1 E€xperiences démands which may be closely 
approximated with a Poisson distribution with 
rate L;. 

Miremarecoe Tepair facility completes repairs on 
items in a constant specified time not greater 
than the procurement lead time t (when repair 
is possible). 


TWiewecocucement: eed time for the depot is a 
G@metbane Value tt. 


Bae -sOOMMOL Cxperlvence crv@ssed orders, i.e. 
emders are received from the depot in the same 
Sequcnee in waileh they were placed. 


The depot is never required to place an order 
come haghtersrevyel Of supply when a previous 
order has not been received. This assumption 
does not usually distort the model since for 
mmc usvel case where ordering costs are high 
and holding costs relatively low, orders are 
placed intrequently. 


An order quantity of one is appropriate to a 
minimum expected cost solution for base i. As 
argued by Feeney and Sherbrooke [1] , where 
holding costs at the lower level are typically 
hagner than at the higher level because of 
fwem tactors as restricted storage capacity 





and ordering costs are much lower at the lower 

teveleeleameat the higher level since they only 

entail administrative and handling costs 

instead of for example, contract negotiation 

costs; a policy which orders small amounts 

often and thus maintains low stock levels 

usually compares very favorably with the mini- 

Mum cegeected cost policy. This is especially 

true where integer quantities only are con- 

sidered. 
B. PARAMETERS AND VARIABLES 

A few comments in explanation are appropriate for some 
of the parameters and variables listed in Table If. Lie ey 
wee teeanad ft, are inventory holding costs for a base i 
stock, base i repair, depot stock and depot repair respec- 
tively. These figures represent the costs composed of such 
things as storage, handling, loss, breakage, etc. for 
faerie weach Lkem OL inventory for a gqiven length of time. 
ft is assumed that the holding costs are expressed in units 
Of dollars per item-year. Therefore, the instantaneous 
Yate of incurring carrying charges can be expressed as Ix 
where x is the on hand inventory level. Although diffi- 
mbt CO Specify in practice and often not varying linearly 
With stock on hand, as assumed here, it is considered 
necessary and acceptable to introduce holding costs in this 
manner. 
The Order cost for the depot, A, is taken to be total 

COst assocaated with placing, following up and receiving an 


order. Components of A might be the costs of letting a 


Contract and receiving the order at the warehouse. 


LO 
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Ly mean rate of demand for the item experienced at 
base 1 


Pi orobability that a carcass turned in at base 1 
may be repaired at base 1 


ler expected repair time for items repaired at base i 

oF mean procurement lead time experienced at base i 

ry base i reorder level 

woe inventory holding cost for base i stock 

Tir inventory holding cost for hase i repair 

cy probability that an item forwarded to the depot 
from base i repair can be repaired at the depot 

ee constant repair time for items repaired at depot 

ie procurement lead time experienced at depot - a 
constant 

16 Gepet reorder level 

Q depot reorder quantity 

I inventory holding cost for the depot stock 

Ir inventory holding cost for the depot repair 

A order cost for the depot 


Bee(ty, 04) probability that base iis able to fill an order 
without delay whenever a demand is experienced 


ayi expected order/shipping c<ime for base i given 
that sufficient stock is on hand at depot for an 
immediate fill - major components of dj are 


shipping time and time required for processing 
of papers, packaging POmmenesoLdeck 


Table I. Parameters and Variables 


MZ 





The SOUUTION OF s2HE MODEL 


Solution of an inventory model for an optimal policy 
is often a very difficult task. First, it is necessary to 
determine the structure of the problem and what information 
is available and desired. Usually the analyst proceeds as 
though the consequences or dollar costs of certain states 
of the system modeled can be determined; examples are the 
Costs of holding an item in stock without issuing it fora 
Wear and the cost of placing an order for some specified 
number of articles. When all such states cen be assigned a 
cost the investigator then tries tc determine the proba- 
bility of the occurrence of each state of the system and 
then the expected cost of maintaining the system fcr 4 
given period of time. It is necessary that this expected 
cost expression be a function of variables which may De 
determined by an inventory policy; otherwise, expected 
operating costs are independent of inventory policy and 
there is no problem. Most commonly, these policy variables 
are lead times, order quantities ang reorder levels. It 1s 
then necessary to solve the expected cost expression for a 
policy which yields the minimum expected cost. Of course, 
expected costs may not be of great importance and dite Weve pe 
be of value to design a policy which, for example, guaran 
tees with a certain probability that all demands received 
can be filled within a given period of time. The analogy 


holds in this case also. The consequences of each state 


3 





of the system in terms of the probability of issue within 
the given period of time must be determined and the probha- 
bility of filling an order within the given period of time 
expressed as a function of variables which the inventory 
policy controls. As in the case with expected costs Lee Ls 
then necessary to solve this problem for a "best" policy 
which in this case is the one which produces the desired 
probability that issues are made within the given time 
period. It is important to realize that the mathematical 
problem which results from the type of formulation discussed 
above is unconstrained; that is, any policy which produces 
the desired value of the Single criteria is acceptable. 

A problem occasionally arises, however, where it is not 
possible to assign casts to each possible state WilLCchware 
Compatible with costs associated with the characteristics 
desired. For example, while an investigator may attempt to 
find policies which minimize total expected system dollar 
costs over some period of time to operate a system, Tel abs 
not always possible to assign costs to factors such aS not 
being able to fill a customer demand. What, for instance, 
is the cost in dollars of not having an item in Seock CO 
repair a ship damaged at sea or a meteor struck space Cla lime 
In these situations the best that the analyst can usually 
hope for is to "pass the buck" as far as assigning costs of 
being out of stock igs concerned to the authority responsi- 
ble for accomplishing the mission involved, i.e. the eel 


sion maker. In this formulation the analyst —-oduces for 


14 





the decision maker a series of alternatives of expected 
costs against probability of being able to fill a demand. 
The alternatives are usually arrived at by solving a group 
of problems where expected value of the sum of the compar- 
able costs is minimized subject to constraints on the non- 
comparable costs or consequences; in this case the proba- 
iombty Of not being able to fill a demand. It is then up 
to the decision maker to select a policy which corresponds 
to his opinions or feelings about the trade off between 
total system expected cost and payoff, or in the parlance 
of systems analysis he must specify his preference ordering 
on the set of measures of cost and effectiveness. If the 
analyst is able to deduce or is supplied with an expression 


for the decision maker's trade off policy between ‘elelone 


no 
4 4% 


~) 


effectiveness in this situation, he may, of course, apply 
Eas policy towehe list of alternatives he has developed 
and determine a "best" or "optimal" policy. It can be 
shown | Hadley and Whitin, 1963 | that in picking a policy or 
formulating a preference ordering the decision maker TS) 
effect placing a dollar value cost on not being able to 
fill a demand and in this sense the two formulations 
discussed above may be considered equivalent. 

In this thesis a formulation similar to the non- 
Comparable costs example discussed above is used. First it 
is necessary to specify a set of costs to be considered; or 
stated differently, it 1s necessary to determine those 


characteristics of a state of the system which may be 


iS 





aaemeea COstSs which are Comparable. Other characteristics 
Maweh arc important but which are not included above are 
then treated as constraining factors and labeled effective- 
mess Criteria. A set of points in the cost-effectiveness 
weeecelnaerced above is determined by the solution of a 
program: 

Meme Zze total expected system costs 

(comparable) subject to constraints which 

requame that each measure of effectiveness 

(non-comparable cost) meets or exceeds 

some stated value 
tem many diffesent levels of the constraints. Each of these 
points is simply an alternative presented to the decision 
maker as in the example noted above. It is then necessary 
to apply the decision maker's preference order, if avail- 
able, and solve the related problem for an "optimal" policy. 

Mitoece (ate ely elaborcatesoptbimization scheme is 

fee oped min oreer tO Circumvenmt the very difficult task of 
meet ty Ingres Cost Of not being able to £111 an order 
Perec lacteus, tie backorder cost. it is an attempt to place 
Pes IneLOoaduction Of this consideration into the system model 
Secetily through the decision maker, via his preference 
Grdering. Therefore measures of effectiveness must neces- 
Sarily address the manner in which backorders are to be 
eaeaeea and more particularly the rate of backorders experi- 


enced. A measure of effectiveness is: 


Bence; , rc.) = themea@odorlity that base 1 is 


able to fill a demand immediately from stock 
on hand. This has been referred to by 
several authors as the fill rate. 


Tiers one such measure of effectiveness for each base. 
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The total expected system cost to be minimized is 
considered to consist of the following comparable costs: 
Perel expected @rder cost for the depot and total expected 
holding costs summed over bases anc the depot. Although 
these components of total system cost are also difficult to 
assess in certain situations, they do not present the same 
possible elements of "operational imperative" as the back- 
Or~aer State as the examples of the damaged ship and space 
Cratt illustrate. For this reason and because of the diffi- 
Mimeweot abkbiving at an Optimal policy via a preference 
ordering or decision maker's judgment when the set of 
measures of effectiveness is large (the decision maker is 
offered too many alternatives) these costs are assumed 
Peo eleodeoncwineliucded in the objective equation of total 
system cost. Mathematical expressions for the components 


Simume program are developed below. 


A. CONSTRAINTS 

Let B,;(t;,r;) be the probability that base i is able 
Peer lit a given demand immediately from stock on hand. 
Pine uu; to be the number of items in repair at base i 
Pam tnOse Items on Order at base i, Vv; to be the number 
of items in repair at base feos te and ws, to be the 
number of items on order at base i. Then, since the inven- 
tory position of each base is a constant r- + 1 under the 


sil 


one for one ordering policy used by the bases: 
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Bi(€,,r,) = P(quantity in base i stock > 0) 


ry 
-» p(uy = 3) 

j= 

rij 
2 gy p( vy = q, wy = j - @) 

fo c=0 

ry j 
= > {ewe piwy= j - a}. 
j=0 g=0 


This result follows from a finding presented by Richards (3 ] 
Peech Shows @hat the distributions of the number of units 

in repair and the nuyser cif units on 
A summary of his results follows. Item failures or aemands 
placed upon the system in an interval of time (u, uts| 
constitute a Poisson process {r(s),s = o} . When a failure 
wears the item returned te base i is designated bee 
repairable with probability p; and is, therefore, not 

base i repairable with probability 1-p;- It is assumed 
that the classification of failed items as base repairable 
Or not base repairable is independent of the process gener- 
ating failures and is also independent among failed items. 
In addition, it is easily shown that the number of items 
designated base i repairable in some interval of ete 

(u,uts] , i ToePOlscOn ween meGan rere Lip, and the 


number of items which are not base repairable in that length 


ie 





of time is Poisson with mean rate lL j(1-p;). Also, the 
(eee ove von Otethe number of items base 1 repairable ina 
length of time given that n total demands have been made 
in the same length of time is Simply binomial with para- 
Pease neecaad p;,. With these facts in mind it is only 
necessary to show that the number of items base i repair- 
able received in a period of time is independent of the 
number of items which are not base i repairable received in 
that same period of time. If this is true, it is well known 
Mite thounumboer of items in each classification received in 


different, possibly overlapping, time intervals is also 


°) 
oa 


independent. Consider the following: 


II 


p { P, (uts) pee ier 6(urSs) = 5 ey, 


II 
lI 


p {Pp (ats) _ F5 (u) =k F (uts) - F(u) 


Pp { P(uts) - F(u) = xem | 


“Lys kt 
eo ( Se Se) Gis) ee 
k 1- pi) Pi (k +m): 


Maa aS (l-p.)L.s 
—  [LO-pi) tis] (pghis) e ~~ e on. 
¥ iss mgule 
K -(1-p,)L.8 m =-—p,Lis 
[(1-pi) Lis] e bow CLs) Soa 
nn ee OE ES 
k. He 


I! 


p {Fe(uts) - F5(u) = bp [Py (us -- Fy (y) = of 


ee, 





Because the number of items on order at any time varies 

only with the receipt of demands accompanied by items which 
are not base repairable and the procurement lead time, the 
Gmcdepericence of the distributions of v, and w, is esta- 
blished. Because of an important theorem by Palm 14] sels 
Peso known tiat the steady state probability distribution 

of the number of items in repair is Poisson distributed with 


Parameter rt L and the number of items on order is 


iPi 
Poisson distributed with parameter t,L; (1-p,). ft is 
Piteresting to observe that these distributions depend on 


the lead times only through their mean values. Thus 


ae, , 
Beea ec; ) =) i? (ptarre ney p(jraitihy (pa) ) 
-O 
Where a . 
= Or 
DCSE Ay 1): ieee =a 


BeeeeObJECTIVE FUNCTION 
Let K(r,Q,r%,,%o) be the total expected system cost bRONG 


a year. Similarly let a 5 ate and K, (r,Q) peste 


i) 
expected costs for operating base i and the depot respec-— 
eet «Then, K(r,Q.r,,r5) = Ky (x, €)) + toe) at K3(r,Q). 


1. Development of the Depot Expected Cost Expression 


Let 
SD(r,Q) = E(number of items in stock at the depot) 
RD = E(number of items in repair at the depot) 
RSD(r,Q)= E(inventory status) 
OD(r,Q) = E(number of orders per year) 
BD(r,Q) = E(number of backorders at the dep: c) 


20 





where E(X) denotes the expected value of the random 
variable xX, and the inventory status is defined to be the 
number of items in stock and in repair at the depot less 
any backordered items. 

TC Can be shown that: 

eerie) — TN) OD(r,Q) + (1,)RD + (I)SD(r,Q). 

Since the average annual demand experienced by the depot 
hs) + £OrKr L, (1-4) -P}) + Lo (1-q5-p,) = L items the expected 
number of orders placed by the depot in a year is L/Q 
where Q is the size of the order placed. It may also be 
shown using Palm's theorem that the expected number of items 
in depot repair, RD, is simply L(rt). Since the inventory 
status includes the number of items on hand plus those in 


- A) - 71 Fr 
repair less backorders, SD(r, 


ie) 


) = RSD(r,Q) = RD + BN(r,Q). 
mamecdaition 1t can be shown that BD(r,Q) = te ee with 


probability one where P is the probability that the 


Out 
number of items on hand is less than or equal to Zero. 
eet sce eropca sin the next section. Therefore, it is 
CGnly necessary to determine the expected value of the 
inventory status, RSD(r,Q). This, however, is exactly the 
expression developed by Hadley and Whitin 12 fer the 
expected on hand inventory for an inventory model which has 
no repairable feature and experiences demands which are 
distributed Poisson. 

Let 2 denote the depot inventory position. Te canbe 


shown that 2 may assume values rtl,...,r+Q. 2 does not 


equal xr for any finite length of time since when @2@ = rt+l 


Za 





and a demand is received (recall this may be interpreted as 
the case when a demand is received without a repairable 
Carcass) an order is placed immediacely and Z2 = r+Q. It 
is also known (calliner, Morse and Simon, 1959 | that the 
inventory position, Z, is uniformly distributed on the 
values rt+l,....,r+Q, i.e. % takes on values r+j with 
probability 1/Q where j=1,2,...A. 

fet xX be the inventory status. Consider the system 
Sele - t Where t is the procurement lead time. Any 
Order placed after u- t will arrive after time u. Any 
Seeer placed at Or before u- t will arrive not later 
Poe Ue tinen; if 42 = rr} at time u- t the probability 
that there are xX units on hand at u is the probability 
that r+ j- x items were demanded in time t if r+j- 
x 20 and zero if r+ j -x <0. Since demand is Poisson 
distributed with parameter Lt the probability that 
r+ j-x demands are received in a period t is 


-Lt r+ 1-x 
e (Lt) J 


————— a lig mee f= x 7 OC 
Coe eae as 


0 if r+j-e x 0. 
Then 


Q 
p(X = x) =) | PO Z= rt+j) p(Z = +3)} 
j=l 


a2 





Pee] —-8x;: bt) if —-@O8< x <€ r+1 


Il 
Ole 
ae 


Pleat bt) ee eed S rt Q. 


II 
O]r 
naic 


j=x-r 
Mnerefore: 
E+O 
RsD(r,Q) = ) p(x =x). 
Zz 
x= 6 
Pemtmen LOllows that: 
r+Q 
sD(r,Q) = ) — p(X=x) = L(rt) + LP out. 
x =—OO 


2. Development of Base i Expected Cost Expression 
Byeemmeragument Wwiich Gxactly parallels the one 


above for the depot it can be shown that: 


ey = } = 
ee a tj) Tuy ea Le (rt. ) “3 de pete. p(X; x) 


fat 
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This is exactly the expression for Ky (xr, Q=1) where the 
term corresponding to the order cost is omitted since it 
does not depend on the variables ri. © and Q, and the 
demand rate and procurement lead time parameters are those 


epeLovrtate for Base i. 
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C. DEVELOPMENT OF THE EXPECTED BASE i PROCUREMENT LEAD 
mney, 

It is necessary to determine an expression for oh, the 
emoccted Base i procurement lead time, which appears in the 
expressions for Kj (r4,t;) and Enea) pe tie procurement 
lead time for Base i can be decomposed into two parts. The 
first, called the order and shipping time, is that time 
Which is required to notify the depot of the base's need 
meus the actual time required to transport the unit from 
the depot to the base once a unit becomes available. This 
time will be assumed to have expected value d,; and is 
independent of the inventory situation at the depot. The 
P-eoncsparct 1S8ene amount Of time which the depot requires 
Eemobtain a Unit to issue. I£ S denotes the latter tine, 
then it is clear that S will be zero if and cnly if the 
@n hand inventeony at the depot is positive. Thus S isa 
Fandom variable which depends on the on hand inventory at 


depot, the items in repair and the items on order at depot. 


The expected Base i procurement lead time is, therefore: 
t. = (claiee: ECS). 
Define the following terms: 


W 


the random number of items in depot repair 


ae) 
tl 


the random number of items in depot supply 
ready for issue less any backorders 


Y =WetdH 
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Pewee eve miuneer Of peceipts by the depot from 
procurement and repair in any time interval 
(u, uts] of length s 


a ee omOnmiwe hor wlLeCccipis soy che depot from repair 
in any time interval (u,uts| of length s 


Z. = the number of receipts by the depot from 
procurement in any time interval (u,uts} of 
length» s. 


It remains to determine an expression for E(S) asa 
memetion Of the program variables r and Q to completely 
Specify the program statement sought. For 5 2 OF ast aS 


true that 


od 
pis€s) = ) pis $s | = -n) p(H = -h) + p(x > 0) 
h=0 


wnere 
CO 
— 
Pets oh) = 7) P{Y-x-h) p(w = x) 
x=0 
and Se 
4 7 
Pout > ie LH ~ = 
h-<« 


As developed in the preceding section it is known that 


Q 
p (Y=k) =< wi (Gomes i Gane k 2 x41 
j=l 
Q 
=~ i \ p(r+j - k; Lt), rt+l ok Sr+t+Q 
= 


in the steady state. Therefore, 


Za) 





rth Q 


ign Q 
: 1 \- 4 
+ bet al — — e e 
5 p | xrt+j (x-h) ;Lt 
x=rtht+l J=saee 


as 
‘s) bs (et) D | 
ar 
where i LyQ) ie L5Q>- 


Also as shown in the preceding section the random variables 


Peeeance Y ame A2andependent. Therefore, it is possible to 


write 
oO CO oo 
S : ail 
p(H> 0) =) p(Hsh) = »  p(vehtx) p(W=x) 
ig reg! QO 
= yj > ‘ a PiGeas a(tlt)G BE) 
hei Sc) j=l 


r+OQ epi O al 3 Q 
p (><; (rey?) } — >. (2 ». 


h=x+l x=r+1-h j=x-h-x 


p(r+j -(xth);Lt) p(x; (ct) D ) f : 
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It remains, therefore, to develop an expression for P(S £5 
H = —-h) 

Case 1: Let s <t. (Recall that t denotes the depot 

procurement lead time.) Whenever H= -h, S_ will be less 


man Or equal te s if and only if R, is greater than h. 


Thus , 
co 
pis4s |H=-h) =p(Rg>h) =) —p(Rg= 5). 
j=htl 


Recall that a necessary consequence of the assumption that 
mememore than one order can be unfilled at any time is that 
only one order can be received in any time interval of 


Yength s<t. ‘Thus, 


p(Ro=3? a p(Z.=0, X.=]3) + p(Z =Q; Xe a apa) 


oo 
o> 


Semcinuing, 
p(Z, = i, X, = j) = p(@s = i) p(X, = §) 


pace the random variables 2, and xX, are independent. 


Proceeding, 


i 4 =O ) = P(no order arrives in an interval (u, uts| Gis 
length § s) 


P(no items are on order at time u+s-t) 


lI 


Tat 


p(y = k) 
ke 41. 
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EO Q 

yy 3 Gey aa; Lt) and 
p abs 

k=r+1 j=k-r 


O]re 


| 
oe 
I 


p(Z.=Q) p( 4, = 0). 


Smmeease 1: Let s<¢ rt. It is true that: 


C~ 
p(x, = 5) = ) p(X = 5 | W= x) p(wex), 
X=] 
Then 
P(x =j | W=x) = P(x-j items are in repair at time uts 


each of which was in repair at time u)_ 


mmeexXpression for the steady state probability that x 
items are in repair each of which has been in repair s 
units of time can be developed for any repair time distri- 
bution as follows [ Takacs 1962 | Se Wiieee 
a(s) = P(any item in repair in the steady state has 
been Iererairedeeteadet S Units of tame): 
Mime from the distribution of age for a renewal process it 


is known that 


CO 
ra 6S) a t / fi - P(t)| Ct 
| S 


A 
where T is the mean repair time and F(t) is the 
peebabrigty distribution function for the repair time. 

Let Dy be the event “"y units are in repair 1t time uts 


Gach Of which was in repair at time wu" and let xX be the 


number of items in repair at time uu. Then 
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Go 


P(Dy) = ) P(Dy | xX=x) P(x=x) 
x= j 
rene) Pa ee 4 a 
= -Ls 
= ) | a(s)¥ (1-a(s))* ¥ 408) ¢ 
=y : re 
A - CO A 
_ a(s)¥ (Is)¥ e % [Zs (1-a(s)) | *-¥ 
y- xy (x-y). 
= 
A —-Ls 
y A 
2 a(s)¥ (Ls) = als (1-a(s)) 
y: 
x “NN 
(als) tel” vomne a(s) A 
a Oe ee ig p(y;Lsa(s) ) 
ivet, 
: ; A 
meee j) | N= x) = Pizeve tisa(s))- 


POretene Case where the repair time is constant, the cumula- 


Mee aistribution function as given by: 


IV 


ii. a ee 


I 
rt 


F(u) 


= 6 Se) we sees 
figererore a(s) = 1 - s/rt. 


Suess 2: Let rts < t. 

Since rt€s everything in repair at time u will have 
Peeg@epeeaired by time wuts. in addition, any items which 
Syeer repair in the interval (u, uts-rt] will also have been 
repaired. Therefore, the conditional probability ese, sclal= 


number of items received from repair in the interval 
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(u, urs] is j given that x items are in repair at time 
u is 
x 


p(X,=j | Wex) = p(j-x; (s-rt) L) if j na te 


For the case in which j < x. 


p'(X,=j | W=x) = 0. 


Therefore, for O€£s < rt: 
OO r+Q Q ; 
p (S€s H=-h) = a "3 a yy p (r+i-k; Lt) |- 
j=ht+l k=r+1 i=k-r 
0° x 
; aN 
» Ie | x3: eeMGEN) | p E (ict ) D | “ 
X= Jj 
a r+Q @) - exe) 
| 4 wr a ‘e PF f+. 3 iv ~ &. | . 
Ha as 7 o Zz bP Nde t ae ak = oe 
k=r+1 I= Ke —T X= j -O 


p [x-j+0; usa(s)] p[ x: (rt) T] 


Ana for rt © sxt: 


oO ( 1a h®, 


Q 
p(sés | H=-h) = on | > 2 y p(x+i-k; bt) | . 
k 


j=nhtl 


j 
> p | j-x: (S-rt) tT | p | 23 (ee) oo + 


0 
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O jk 


6419) Q =O 
E - ) )  -p(rti-k; 2 | yy 


k=r+1 i=k-r i x=O0 


A) 
L 


p | 5-0-x: (s-rt) L | p | x: (at) t | 


Meace 2: Let sa =e 
Because it is assumed that not more than one order is 
outstanding at any one time it follows that: 
p (Ss | H=-h) = 1. 
With the expression for P(S£s) thus developed it is 


Possible to calculate: 


E(s) = f p (s> s)ds . 


D. STATEMENT OF THE PROGRAM 
The minimum expected cost solution for given values o£ 
Poeeeenia b> where b, is the probability that an order 


Seeurring at Base ican be filled immediately may be deter=- 


termed by solving the nonlinear program: 


Z 
Minimize Ry Ke (xr. t.) + K, (xr,Q) 
ot ee eat d : 
i=l 
subject to 
B,(t,, x) Mae een cy 


where r,.r and Q may only assume non-negative 


mmeecqer values. 
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IV. SOLUTION OF THE PROGRAM 


This thesis does not develop an efficient procedure with 
which the program specified can be solved. However, several 
observations are offered below in an attempt to delineate 


Ene scope and nature of the difficulties involved. 


A. PROGRAM TYPE 

It should be noted that the program developed is 
Hemraneaer in both the objective function and the constraints 
and the program independent variables are ye U5, © and Q, 
eacn of which may assume non-negative, integer values. It 
is not apparent that the program is convex, and it is almost 
Seiseadierrnat establishing that the program is or vis not 
convex in the objective equation and/or constraints would 


involve laborious computations. 


B. DEVELOPING A SOLUTION PROCEDURE 

ie is important to note that for given values of YF and 
O, ©; and ry are determinate, i.e. they assume the smallest 
integer values such that Belen a Gey) > b, and Boers) op 
because the objective function which is to be minimized is a 
non-decreasing function of fr, and rn as are the values 
of the constraint expression. This indicates that an opti- 
mum solution can be arrived at by varying r and Q alone 
and allowing ry and E>) eO@sussumewene Smallest ineege: 


Values such that the constraint expressions are met. 
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It 1s reasonable to suppose that due to the nature of 
the computations involved in the solution of the program, 
for example, summing various expressions which are easily 
located in tables, a high speed, digital computer would be 
employed. While knowledge of convexity can be extremely 
metptul in locating optimal solutions and in fact assures 
Ewe any local optimum is global for the problem, the 
limited nature of the search region and certain special 
features of the program give grounds for optimism that the 
Ppmegiam can be solved rapidly. 

It should be noted that the bulk of the computations 
involved in a single iteration of a computer search routine 
involve determination of the t, 8 as EUNCtzOnS of © famed 
Oeics Probably the case that many of the terms in the 
expression for en could be shown to be so small relative 
to the others that they could be disregarded without loss of 
essential accuracy. Central to the computation of t. 1s 


the expression 


| ‘@) 
s 


P(S€s) = P(S€s | H=-h) P(H=-h). 
h=0 


While at first glance this appears to involve formidable 
computations, the physical interpretation of H_ should not 
Peer emliooked. dtvis highly doubtful that any optimum solu- 
tion to the program would be characterized by values of fr 


and Q which would produce a significant probability that 


H be smaller than, for example, -3 or -4, depending on the 
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program parameters involved. Certainly large probabilities 
of small values of H would increase t, to a point where 
Mmemvalues of r), and C5 See reds ee sewer! thercon= 
Straints would be prohibitively large, especially when as is 
usually the case the holding costs at the bases are larger 
than those at depot. In this example it follows that compu- 
metwon of P(S£s) toa degree Of acCuUbacy SHEETcCIent for 
the purposes of the program would involve summing over a 
very few values of the dummy variable h since the remain- 
ing terms would be extremely small. 

Consider also the regions of r and Q which should be 
eeerched. There is reason to suppose that the ranges of 
femes Of r and QO which produce minimum expected cost 
Zemmerous tO the program are restricted. Certainly for some 
value of r (even if Q=1) the contribution made by adding an 
Peemetonal item to the reorder point does not significantly 


decrease the values of the. t,'s and, therefore, larger 
Values of r need not be considered. A similar argument 
can be made for Q. There is a value of Q such that the 
benefit in decreased ordering cost gained by increasing the 
value of Q does not offset the penalty incurred by 
increased holding costs which also result. These border 
values on the regions which are searched for optimum solu- 


mrons to the program are, Of course, extremely sensitive 


tO tite program paramewer values. 
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V. APPLICATIONS OF THE MODEL 


As neues when developing the model it is often the 
objective of a constrained optimization problem to develop a 
Group of alternatives from which a decision maker can select 
a "best" element which reflects his feelings about the trade 
off between comparable and non-comparable costs. In this 
thesis an analysis has produced a program which provides a 
method to obtain points in a space of costs and measures of 
effectiveness. For a particular inventory item and situa- 
Cion such a set of points might be as shown in Fig. 2 where 
b; xrepresents the probability that Base i is able to fill 
an order on demand, and associated with each set of values 
of by and by is a cost rounced to the nearest thousand 


eo llars. 


i ee eel eli Ea ia 1S 15 sl 
=| is 
ie) 


1 ae a 
6 Cee oe 13 
9 

bo 


Figure 2. A Representative Set of Decision Alternatives 
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While it is not the goal of this thesis to investigate the 
various methods which can be used to develop a cost- 
effectiveness surface, the following points are worthy of 
Mo ce : 

The logic of the» physical realities of the model 
dictates that costs will be non-decreasing in Db) Semd “bo. 
If this were not the case it would be possible to “buy" a 
greater amount of effectiveness for less than is needed to 
pay for some lesser degree of protection which is contradic- 
Eoay to the minimum cost result of the efficient solution. 

In many situations it will only be necessary to develop 
minimum cost solutions for small ranges of possible values 
bb) and b>. The reason fer thc 1s thae  fimemeny 
instances the decision maker can specify in advance, without 
immediate knowledge of the consequences in terms of effec- 
tiveness, the maximum expected cost which can be incurred 
and/or the general range of effectiveness he wishes to 
echaeve. It may be noted that certain ranges of effective— 
ness might be of particular interest and that for those 
areas the effect of a smaller change in the valves of by 
and by could be sought; therefore, it would be necessary 
Femcearch along a finer grid of walues for 6, jand Doeeee 
tease for some particular segment,_of their ranges. Ie Can 


be said, therefore, that the ranges of b and by», to be 


1. 
searched can be narrowed considerably with knowledge of the 
decision maker's general policy and the physical realities 


Siena cystem to which the model 1s "ap plied: 
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It might have been noted that terms in the objective 
equation involving the repair facilities' holding costs 
were independent of the program variables and, therefore, 
Could have been deleted. These terms were included, how- 
ever, in order to show the ease with which this model and 
Peegram can be employed to discover not only efficient 
inventory ordering policies but also efficient inventory 
Sisaering and repair facility management policies. Suppose, 
for example, it is possible, by incurring an extra cost, to 
decrease the repair time at the depot by some fraction of 
the original repair time, that is there exists a function 


C(rt) defined for rt in some interval [a,b], which is 


decreasing in rt. It is then possible to specify a new 
program: 
Minimize ss ae ee O, Yvert Cc’ Gse) 
subject to B; (ty, x.) > b, 
vim 826 a 
=—rt |= = 
i= 1, 2 


where Y), Yor © and Q are non-negative 


integers and rt is some real number. 


Similar cost terms and comereatnte can be added for base 
repairs or any combination of repairs where the repair times 
can be considered variables. In the case of the bases, 
however, the additional costs need only be considered to 
vary with the expected base repair times. It is then pos- 


sible as above to present the decision maxer vw.th a revised 
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set of costs corresponding to the various degrees of effec- 
tiveness achieved. 

It may also be noted that while the model configuration 
developed has two bases, there is no reason why systems with 
any number of lower echelon units cannot be modeled. The 
Only changes required in the program formulation would be 
to increase the number of cost terms in the objective 
function and the number of constraints by a number of terms 
equal to the additional bases modeled. It may be observed 
Enact the model specified can be weed to deseribe a multi- 
echelon system where item repair is not considered. This is 
possible by stipulating the probabilities of repair at the 


Various facilities to be zero. 
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